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A  study  of  Dynamic  Systems  with  three  variable  para- 
meters is  made  by  an  initial  scanning  of  the  basic  geo- 
metric properties  on  the  three  dimensional  space  generated 
by  these  parameters.   From  these  geometric  properties,  a 
Root-locus  technique  that  simplifies  greatly  the  amount  of 
work  in  analysis  and  design  is  developed.   This  technique 
is  extended  to  systems  with  k  variable  parameters. 
Finally,  singular  lines  in  the  parameter  plane  are  treated 
as  a  special  case  of  the  parameter  space,  and  formulae  are 
derived  for  a  fourth  order  system  leaving  the  field  open 
for  systems  of  higher  order. 
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1.   INTRODUCTION  OF  THE  THIRD  PARAMETER 

This  study  takes  as  a  starting  point  Mitrovic's  ap- 
proach to  the  root-finding  problem.   Although  this  is 
essentially  a  two-variable  method  which  will  yield  only 
two  relations  out  of  the  three  needed  for  the  three-para- 
meter case,  the  parametric  equations  obtained  will  help  in 
understanding  the  nature  of  the  basic  geometric  properties 
of  the  parameter  space  generated  by  &,    j9,  and  y,    the  three 
variable  parameters.   The  fact  that  it  is  possible  to 
obtain  only  two  relations  is  justified  because  a  point  on 
the  s-plane  is  completely  defined  by  only  two  conditions, 

and  therefore  any  other  possible  relations  will  be  linear- 

2 
ly  related  to  the  two  originated  by  Mitrovic's   root-find- 
ing method. 

Let  us  now  describe  briefly  Mitrovic's  approach,  as  a 
base  for  the  present  study: 

The  characteristic  equation  for  a  given  system  can  be 
written  in  a  general  form  as: 

F(S)  =  a  Sn  +  a   .  Sn-1  + +  a^S  +  a   =0        1.1 

n      n-1  1     o 

or  in  a  closed  form  as: 


n 


F(S)  =  \   avSk 


k=0 
If  now  S  is  expressed  as: 


S  =  -Ceo  i  jU)Jl-Z2  1.3 


n    -"    n 


where  C  and  CO  are  defined  in  figure  1.1,  which  is  the  re- 
presentation of  a  point  in  the  s-plane  that  is  defined  by- 
assigning  values  to  C  and  u)   . 


0  = 


cos   (-C) 


Fig.  1.1   -   Representation  of  a  point 
on  the  s-plane 


Equation  1.3,  written  in  another  form  yields: 
S  =  O)(cos  0  +   j  sin  0)  =  05e- 


J  9 


1.4 


—  1  k 

where    0   =   cos       (-C)  .      Then      S      can  be  written: 

Sk   =   dV1*0   =  wk(cos   k0   +   j    sin  k0)  , 


1.5 


from  which  the  Chebyshev  functions  are  obvious  and  are 
T,  (C)  =  cos  k0  =  cos(kcos   C) 


1.6 


77    /f\         sin  k0    ..   sin(kcos      C) 
k  sin   0  sin(cos~*    £) 


1.7 
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S  =   us 


(-C)   +  j  Vl^Fu  (-£)] 


=  o)k[(-i)k  T]c(C)  +  j  7i^P(-i)k+1uk(C)] 


1.8 


1.9 


Inserting  these  in  1.2  and  requiring  that  the  reals  and 
imaginaries  go  to  zero  independently,  provides  the  two 
equations : 
n 


I  aka,n*(-l)kVO  =  0 


k=0 


n 


1.10 


k,  ,%k+l. 


Zvn  f-1'     V?>  =0 

k=0 


But  t^o  =  ^(f3-(Vi(s*) 


and  upon  substitution,  the  following  equations  are  obtained 
n 


Ic-dVJi*  uk-i(C)  =  ° 


k=0 
n 


1.11 


li-l)\u£   uk(C)  =  o 
k=0 


where  U,  can  be  obtained  from  the  recurrence  formula: 
k 


where  U  (C)  =0,  and   U-^C)  =  1 

Now,  the  coefficients  are  defined  to  be: 


ak  =  Qbk  +  ^ck  +  ydk  +  fk 


1.12 


1.13 
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where  a.,    0,  and  y  are  the  independent  variable  parameters 
of  the  system  under  analysis.   By  substituting  this  ex- 
pression for  a,  in  equations  1.11,  these  equations  can  be 
rewritten  as  follows j 


1.14 


ttB,  +  £c   +  yD,  +  F,  =  0 

olb2   +   Pc2   +   yD2  +  F2  =  0 

where 

n  n 

Bl  =  k-»\   <"nVl<C>    B2  =  I{-»\   <°nkVC> 

k=0  k=0 

n  n 

ci  =  h-»\  »„Vi(»    c2  -  h-v\  %\™ 

k=0  k=0 

1.15 
n  n 

Dl   =     k-»K    "n\-l<?>         D2   =     &-l>\   WnVC> 

k=0  k=0 

n  n 

k=0  k=0 

Information  concerning  the  design  properties  of  the  system 
must  be  obtained  from  these  expressions  (equations  1.14). 
This  information  can  be: 

a) .   For  a  given  set  of  Performance  Specifications, 
what  values  of  ol,    j3,  and  y,  the  variable  parameters  of  the 
system,  are  needed  in  order  to  obtain  the  dominance  con- 
ditions that  will  allow  us  to  meet  the  given  performance 
specifications? 
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b)  .   The  physical  meaning  of  the  values  of  a.,    &,    and 
y,  whether  or  not  these  values  are  actually  attainable, 
and  if  not,  what  modifications  must  be  introduced  to  the 
system,  such  that  it  is  possible  to  obtain  the  desired 
values. 

c) .   Existence  of  dominance  conditions,  if  these 
exist,  and  the  locations  of  the  remaining  roots  on  the 
s-plane,  and  their  influence  in  the  behavior  of  the  system. 

If  C  and  Ui     are  given,  then  the  set  of  equations  1014, 
will  contain  three  unknowns :  OL,    j9,  and  y.   Therefore  the 
surfaces  generated  on  the  a,  £,  y  coordinate  system  can  be 
expressed  only  in  a  parametric  form,  generating  sets  of 
parameter-planes  in  a  and  j8  for  a  given  value  of  y;  or  it 
is  possible  to  obtain  a  three  dimensional  representation 
which  is  not  a  practical  tool  for  engineering  purposes. 
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2.   REAL  ROOTS 

The  characteristic  equation  for  a  system  can  be  re- 
written as 

OA(S)  +  £b(S)  +  yc(s)  +  D(S)  =  0  2.1 

By  making  S  =  -  (7,,  A(S),  B(S),  C(S)  and  D(S)  ,  become  real 
numbers  A,  B,  C  and  Df  so  that  the  above  equation  becomes 

OA  +  j3B  +  yC  +  D  =  0  2.2 

which  is  the  equation  of  a  plane  on  the  system  of  coordi- 
nates generated  by  QL   ,    ft,    and  y,    with  intercepts  as 
follows : 

on  the  OL   axis  oi  =   -  D/A 

on  the  ft   axis  ft   =  -  D/B  2.3 

on  the  y   axis  y  =   -   D/C 
This  is,  then  the  first  set  of  known  surfaces  that  can  be 
used  in  the  three  dimensional  space  generated  by  OL,    ft,    and 
y .   However, there  is  need  for  a  graphical  representation 
which  allows  the  determination  of  the  desired  information 
without  tedious  use  of  computarized  data,  with  which  it  is 
hard  to  get  a  clear  vision  of  the  trend  of  the  different 
values  for  changes  in  the  parameters  of  the  system.   In 
order  to  obtain  this  graphical  representation,  it  is  possi- 
ble to  use  descriptive  geometry  as  a  useful  tool  for  the 
three-dimensional  analysis. 

The  construction  rules  are  simple  and  easy  to  retain 
and  with  some  practice  it  can  be  a  tool  that  can  be  used 
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very  often  in  this  study. 

The  first  step  is  the  representation  of  a  plane  in  a 
three-dimensional  space:   A  plane  is  completely  defined  if 
two  lines  belonging  to  such  a  plane  are  defined.   It  is 
possible  to  draw  these  two  lines  by  determining  the  inter- 
cepts of  the  plane  which  represents  S  =  -  or,  in  the  ot,    )3, 
y  space,  with  the  a/3  and  Oiy   planes.   (See  figure  2.1a  which 
is  a  representation  of  the  plane  S  —  -  <T.  on  the  first 
octant  of  the  coordinates  axes) .   Now,  by  rotating  the  y 
axis  in  the  plane  j3y  in  the  manner  shown  in  figure  2.1a, 
and  in  amount  such  that  the  three  axes  lie  in  the  plane 
generated  by  the  axes  ot   and  j3,  figure  2.1b  is  obtained, 
which  contains  all  the  information  concerning  the  three- 
dimensional  graph.   However,  this  time  it  is  possible  to 
work  with  real  values  and  to  plot  all  the  information  in 
one  plane.   Of  course,  once  it  is  learned  how  to  work  with 
this  construction  method  it  will  not  be  necessary  to  have 
the  three-dimensional  graphs.   Let  us  assume  we  are  to 
investigate  whatever  values  of  Ot   and  £,  are  obtainable 
with  given  values  of  S  =  -  (7.  and  y   =  y.  .   In  the  three- 
dimensional  picture,  we  would  pass  a  plane  parallel  to 
the  plane  otfi,    through  the  given  value  of  y   —  y,  in  order 
to  obtain  the  plane  PNM.   As  can  be  seen,  this  is  not  a 
very  practical  method  even  if  there  exists  the  possibility 
of  building  an  actual  three-dimensional  representation  of 
the  plane  S  =  -  a, .   By  the  construction  method  proposed, 
what  has  to  be  done  is  to  project  the  chosen  value  of  y 
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LJi 


S=-ff,  PLANE 


'%         <* 


hOl    ' 


Fig.  2.1a 


Fig.  2.1b 
Fig.  2.1a  -  Representation  of  a  plane  in  three  dimensions 

Fig.  2.1b  -  Representation  of  a  plane  using  projections  and 
all  coordinate  axes  in  one  plane 
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onto  the  line  ab,  determine  the  point  N,  and  project  this 
point  N  onto  N'f  and  finally  draw  a  parallel  line  to  ac 
through  N'.   Then  the  line  de  will  be  the  actual  represen- 
tation of  the  line  MN  of  the  three-dimensional  graph,  from 
which  we  can  read  the  values  of  ot   and  £.   This  construction 
method  can  be  more  easily  understood  by  the  use  of  a  simple 
example : 

Example  2.1;   Let  us  assume  that  a  system  has  the 
following  characteristic  equation 

S4+  (30H-4/3)S3+  (6Q!+3j9)s2+  (10j8+50y)S  +  (5QH-5y)  =  0    2.4 
which  can  be  rearranged  in  the  form  of  equation  2.2  as 
follows ; 

G(S)  =  (3S3+6S2+5)a  +  (4S3+3S2+10S)  j8  +  (50S+5)y+S4  =  0 

If  S  =  -5  is  to  be  investigated,  then  for  this  value  of  S; 

G(S)  =  -220tt  -  475j9  -  245y  +  625  =  0 

and; 

625 
the  OL   intercept  =  -   o?n  =  2.84 

/*  or 

the  j3  intercept  =  -  _/,lc.    =  1.315 
and  the  y   "     =  -  •=|||  =2.55 

Figure  2.2  is  the  representation  of  the  positive  axes  ol,    j9, 
and  y.   The  negative  axes  can  also  be  identified  by  label- 
ing the  )3  axis,  as  the  -y  axis,  the  y  axis  as  the  -j9  axis. 
Lines  ab  and  ac  are  the  intercepts  of  the  plane  S  =  -5 
with  the  a/9  and  ay   planes;   dotted  lines  are  extensions  of 
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^  (borC-r) 


V  IT  or(-/5) 
Fig.  2.2  -   S  =  -5  lines  in  the  OL$y   space  for  Example  2.1 
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the  lines  on  the  negative  side  of  the  axes.   Let  us  assume 
that  we  wish  to  investigate  what  values  of  ot   and  j9  are 
suitable  for  a  given  value  of  y,    say  y  =   -5.   By  setting 
y  =   -5  and  projecting  this  value  of  y   onto  line  ac,  point 
d  is  obtained  from  which  a  line  parallel  to  ab  is  drawn. 
This  is  line  de,  and  this  line  contains  the  information 
about  the  values  of  Qt   and  £  that  will  satisfy  equation  2.4 
for  S  =  -5  and  y  —   -5.   As  a  check,  let  us  choose  the  point 
f;   at  this  point  ot  =   +5,  fi  =   +1.58  and  y  =   -5,  which  sat- 
isfies equation  2.4  if  S  =  -5. 

Another  interesting  feature  of  this  construction  me- 
thod is  that  due  to  the  linearity  properties,  once  the 
"advance"  over  the  a  axis  for  a  given  change  in  y   is  deter- 
mined, a  grid  can  be  easily  drawn  without  the  need  of  pro- 
jecting each  value  of  y,  by  just  advancing  equal  distances 
on  the  CL   axis  and  drawing  lines  parallel  to  ab.   The  loci 
for  y  =  -5.0,  y  =   -7.5,    y   =+2.5  and  y   =  5.0  were  drawn  by 
this  method.   Similar  results  can  be  obtained  by  determining 
the  "advance"  in  the  j3  axis.   In  summary,  all  the  necessary 
information  is:   the  intercepts  of  the  plane  with  the  axes, 
the  "advance"  for  a  desired  increment  of  one  of  the  para- 
meters.  Then  the  surface  generated  by  S  =  -5  will  be  com- 
pletely described.   Therefore,  for  a  given  value  of  S  =  -O-ii 
the  graph  will  look  like  the  one  in  figure  2.3.   When 
several  values  of  (J   are  needed,  confusion  can  be  avoided  by 
identifying  each  pair  of  lines  by  the  common  point  on  the 
Ot   axis  (See  figure  2.3)  . 
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Fig.  2.3  - 


Grid  of  loci  for  several  values  of  real  roots 
in  the  Otfiy   space. 
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This  method  of  construction  will  be  useful  in  cases  in 
which  it  is  desired  to  analyze  well  damped  systems  and 
also  as  an  easy  way  to  construct  a  grid  of  values  of  real 
roots,  which  can  be  superimposed  over  a  previously  deter- 
mined set  of  constant  C  curves,  in  order  to  determine  the 
possible  values  of  real  roots  expected  for  given  values  of 
a,    j3,  y,    C  and  to 

Example  2.2:   Consider  a  voltage  regulator. 
The  compensation  schemes  have  been  worked  out  and  this 
example  will  be  analyzed  through  this  study  so  we  can  ob- 
tain a  better  picture  of  the  use  of  this  construction  me- 
thod. 

After  manipulation  of  the  block  diagram  the  character- 
istic equation  of  the  system  is: 

S4  +  (2xl06K1  +1350)  S3  +  (3xl08K10+4xl08Ki;L+3.85xl05)S2 

+  (1010K,n+2xl010K, ,+36xl06)S  +  2xl010K  K  +2.1xl010  =  0 
v     10        11  a  p 

by  letting  K,  0  =  OL,    k,.,  =  P,    K  K   =  y,  and  rearranging  the 

characteristic  equation  in  the  form  of  equation  2  2,  it 

becomes : 

(2xl06S3+3xl08S2+1010S)a  +  (4xl08S2+2xl010S)  i3  + 

2xl010y+S4+1350S3+3.85xl05S2+36xl06S+2.1xl010  =  0 

For  S  =  -1  the  equation  of  the  plane  is : 

-9702a  -  19600)3  +  20,000y  +  20964.38  =  0 
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*.     4.-U   ~  •  4-      4.      20,964.38    0  1£ 
where  the  OL   intercept  =  -  _' =  2.16, 

4.1,   fl  •  4.      4.      20,964.38    ,  -"- 
the  p  intercept  =  -  _  '  6 —  =  1.06, 

20  964  38 
and   the  y  intercept  =  -  20000    =  _1-04- 

For  S  =  -5,  the  equation  of  the  plane  is: 

-4275a   -  90000   +  2000y   +  2082.95  =  0 
the  OL   intercept  =  -  .„' =  0.488, 

the  p*  intercept  =  -   _QQ*    =  0.231, 
and    the  y   intercept  =  -   onnn =  -1.04. 

The  results  of  this  problem  can  be  easily  checked  on 
figure  2.4  for  the  case  let  us  take  say  a  =  0.5,  j3  =  0  .  5 
by  drawing  a  line  parallel  to  ab  through  the  point  (0L   =  0.5, 
j3  =  0.5),  we  get  point  c  on  the  OL   axis,  and  then  by  project- 
ing this  point  onto  the  line  corresponding  to  cr  =  -1.0  on 
the  plane  ay,    we  get  y  =   -0.315  these  values  satisfy  the 
characteristic  equation  of  the  system  for  S  =  -1.0.   Another 
variation  can  be  obtained  by  selecting  a  value  of  y,    this 
is  done  for  cr  =  -5.0  in  which  y  =   1.5  was  selected,  by  pro- 
jecting and  drawing  a  parallel  line  to  de  we  obtain  the 
line  labeled  (cr  =  -5.0,  y  =   1.5).   Any  pair  of  values  for  ot 
and  j3  on  this  line  will  satisfy  the  characteristic  equation 
for  the  system.   We  have  picked  the  point  B,  which  yields 
a  =  2.0  and  £  =  -0.385  which  checks. 
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tybor(-t) 


0.5 -* 


0.5 


Or  or(-fi) 


Fig.    2.4    -      S   =    -5.0    and   S    =    -1.0    loci    for   Example    2.1 
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3.   COMPLEX  CONJUGATE  ROOTS 

Usually  dominance  of  a  certain  pair  of  complex  roots 
is  desired  in  the  system.   By  dominance  it  is  understood 
that  the  system  will  behave  in  such  a  manner  that  its  dy- 
namic performance  can  be  approximated  by  a  second  order 
system;   that  is,  that  the  influence  of  the  remaining 
roots  of  the  characteristic  equation  will  be  very  small 
on  the  transient  behavior  as  compared  with  the  influence 
of  the  pair  of  dominant  poles.   This  feature  will  allow 
the  designer  to  obtain  a  certain  degree  of  control  towards 
the  achievement  of  some  performance  specifications  such  as 
settling  time,  natural  frequency  and  damping  coefficient 
(overshoot  or  undershoot) .   For  this  analysis  equations 
1.11  can  be  rewritten  as: 

n  n 

k=0  k=0 

n  k 

+  7  1(-Dk\*nkuk-1^  +  Z^VnVl^  =° 
k=0  k=0 

3.1 
n  n  n 

k=0  k=0  k=0 

n 

+  ll-»\«£\M   =  ° 

k=0 
where  a,  has  been  replaced  by 

ak  =  ^k  +  ^Ck  +  rdk  +  ek*  3*2 
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By  assuming  given  values  of  C  and  cc     for  a  desired  domi- 
nance condition,  equations  3.1  can  be  written  as: 
n  n  n        n 

01  ZbkDk  +  *  l?A  +  y  }?A  +  KDk  =  ° 

k=0  k=0  k=0  ]-.=0 

n  n  n  n  3.3 

01  ZbkEk +  *  IckEk +  y    IdkEk +  KEk =  ° 

k=0        k=0         k=0      k=0 


where    Dk  =  (-l)kttnk  Uk-1(C) 


and     R  =  (-1  ku>  k  U,  (C) 
k         n   k 


3.4 


By  evaluating  the  quantities  under  the  summation  sign  these 
two  equations  can  be  further  reduced  to: 


ttB,  +  0C,  +  y-Q      +  E,  =  0,  and 


0tB2    +    Pc2    +   yD2    +   E2  =  0. 


3.5 


These  are  the  equations  of  two  planes  in  the  space  0L,     /?,  y. 

The  locus  of  the  values  of  OL,    fi   and  y   which  satisfy  the 

desired  values  of  C  and  id     is  then  the  straight  line  defin- 

n  ^ 

ing  the  intercept  of  the  two  planes.   It  is  known  that  a 
real  root  can  be  represented  as  a  plane  in  the  space  OL   $  y, 
so  now  we  can  determine  the  intercept  of  the  real-root 
plane  with  the  straight  line  defined  by  the  intercept  of 
the  planes  described  by  equations  3o5.   There  are  two 
possibilities:   a)  the  line  is  not  contained  on  the  plane 
and  then  the  intersection  will  be  a  point  in  the  space  cl   j9 
y;      or  b)  the  line  is  contained  on  the  plane  and  then  we 
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have  a  set  of  values  of  d   j3  y     on  the  line  that  will  ^ 
satisfy  the  desired  conditions.   In  other  words,  it  will 
be  possible  to  solve  a  set  of  three  equations  with  three 
unknowns  as  follows: 
From  equation  3.5 


OB     +   j8c,    +   yD,    +  E1   =  0 


ol&2   +   0C2   +   yD2   +  E2   =  0 


and  from  equation  2.2 

OA  +  0B  +  yC  +  D  =  O 

By  Cramer ' s  rule : 

-£^^^(020-602)  +  E2(C1C-BD1)  -  D(C1D2-C2D1) 

«  = s 


3.6 


B1(-E2C+DD2)  -  B2(-E1C+D  Dj  +  A  (-E1D2+E2D;L)       3  _ 
P   =   a > 


and 

y  = 


B1(-C2D+B  E2)  -  B2  (-C-j^D+BEj^)  +  A  (-C1E2+C2E1) 


where   A=  6^^(020-602)  -  B2(C1C-BD;L)  +  A(C1D2-C2D1) 

Let  us  now  illustrate  this  use  of  known  surfaces  in  the 
OL   j8  y      space  by  an  example : 

Assume  a  third-order  system  with  the  following 
characteristic  equations 

S3  +  (3a+4#)S2  +  (4Q!+3j9)s  +  20!  +  y  =   0  3.8 

Assume  now  that  we  wish  a  pair  of  complex  poles  with 
second-order  dominance  at  S  =  -5  -  8.66j  and  a  third  pole 
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at  S  =  -50,   What  values  of  a,  j3  and  y   must  we  use  to 

obtain  the  desired  pole  values? 

By  applying   S  =  -50.0  to  equation  3 . 8  we  obtain: 

7,302a  +  10,0000  -  149y  -  125,000  =  0       3.9 

and  by  substituting   S  =  -5  -  8,66j  in  equation  3 . 8  we 

obtain : 

298a  +  400j8  -  v  -  1,000  =  0 

3.10 
and       260a  +  400j8  -  30y  =  0 

By  solving  the  system  of  equations  formed  by  3.9  and  3.10 
we  get  values  of 

a  =  7200,  )8  =  -5.385   and  y   =  -9,400 

By  substituting  these  values  in  equation  3.8,  we  get 

S3  +  60S2  +  600  +  500.  =  0 
or         (S  +  50) (S2  +  10S  +  100)  =  0  3.11 

or         (S  +  5;)  (S  +  5  -  8.66j)  (S  +  5  +  8.66J)  =  0 

We  see  that  we  have  three  degrees  of  freedom  and  therefore 
have  the  opportunity  of  choosing  three  of  the  roots  of  the 
characteristic  equation  for  systems  of  order  higher  than 
three,  once  we  have  made  our  choice  of  three  of  the  roots. 
The  remaining  roots  will  depend  on  the  values  of  »,  $  and 
y  (except  in  the  very  special  case  of  situatim  ) previously 
described,  in  which  we  can  find  another  root  by  determining 
the  intersection  of  the  line  with' the  plane  generated  by 
this  extra  root) .   Therefore,  we  have  the  freedom  of  choos- 
ing three  of  the  roots  of  the  polynomial.   However,  for 
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polynomials  of  a  degree  higher  than  three,  this  method  does 
not  tell  us  anything  about  the  remaining  roots  of  the  poly- 
nomial, and  once  the  choice  has  been  made,  the  remaining 
roots  will  be  determined.   So  far,  as  in  the  case  of  the 
parameter  space,  the  designer  has  the  freedom  of  choosing 
several  of  the  roots  of  the  characteristic  equation,  but 
is  unable  to  determine  whether  or  not  the  remaining  roots 
will  help  in  design,  or  even  make  the  system  stable.   This 
is  a  tedious  trial  and  error  process,  and  it  is  very  likely 
that  for  the  chosen  roots  there  will  not  be  a  satisfactory 
set  of  remaining  roots.   There  are  some  useful  methods  like 
the  Polak  method  and  the  Ross-Warren  method,  employing 
filter  compensation,  but  they  are  essentially  one-variable 
parameter  methods  and  are  not  applicable  in  a  general  sense, 
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4.   CONSTANT  C  SURFACES 

It  was  shown  in  Section  3   that  the  locus  on  the  a.   j9  y 
space  for  given  values  of  C  and  u)     was  the  line  defined  by 
the  intersection  of  the  two  planes  represented  by  equations 
3.5  which  are: 


ttB,  +  )9c   +  yD,  +  E   =  0 
and       «B2  +  £c2  +  yD2  +  E2  =  0 


3.5 


The  above  statement  can  be  interpreted  also  as  "for  a  given 
C,  the  constant  to  loci  on  the  o  |8  y  space  are  straight 
lines",  so  that  by  keeping  C  constant  for  each  value  of  U) 
there  will  be  a  set  of  two  planes  in  the  form  of  equations 
3.5  with  the  intercept  a  straight  line  in  the  three  dimen- 
sional space  &   j?  y.   By  projecting  this  straight  line  onto 
one  of  the  planes  generated  by  the  coordinate   axes  we  will 
be  able  to  summarize  all  the  information  pertaining  to  this 
straight  line  in  one  plane.   Figure  4.1  shows  the  intersec- 
tion of  two  planes  ABC  and  DEJ  as  the  line  GH  and  its 
extension  (dotted  line) .   The  projection  of  GH  on  the  plane 
is  the  line  HF  and  its  extension.   As  can  be  seen,  equal 
increments  in  y  will  produce  equal  advances  along  the  pro- 
jection.  Therefore,  linear  interpolation  can  be  done 
along  the  projection  line  HF,  and  the  information  concern- 
ing the  line  of  intersection  can  be  summarized  in  the  a   /9 
plane.   The  line  HF  and  its  extension  can  easily  be  deter- 
mined by  finding  the  points  H  and  F.   Setting  y  =  0  in 
equations  3.5,  and  solving  for  a  and  £  we  obtain  the  point 


29 


PROTECTION  Cf  IMKRSlCTlON LIME 

OH  flft  PLANE 


i 


'I 


Fig.  4.1  -  Intersection  of  two  planes  in  a  three 
dimensional  space. 
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H,  applying  Cramer's  rule: 

(-EC      +    E   Cl) 

<™ m 


and  /9h  = 


4.1 
(-B1E2    +    B2E1) 

AH 


where  AH  =   ^>-S-2    -    B  C 

Similarly  for  point  F,  setting  P   =  0 

-E  D9  +  E?D, 

aP  =  — if —  4-3 

and       0F  =  0.0 

It  is  possible  also,  to  determine  a  "measuring-stick" 
along  the  line  HF  by  determining  the  value  of  y,    given  by 

(-B,E9  +  B~E,  ) 

v   =  — 4  4 

71        AF  *"* 

where     AF  =  B;,D2  =  B2D,  4.5 

As  can  be  seen,  for  an  increment  Ay  =  y    ,    we  must  advance 
along  the  projection  line  equal  distances  HF.   We  can  draw 
the  projection  line  in  the  &   fi     plane  and  by  using  linear 
interpolation,  label  along  this  line  the  values  of  y. 
Figure  4.2  illustrates  the  construction  method  used  in  the 
two  planes  shown  in  the  three-dimensional  picture  in  Fi- 
gure 4.1.   For  this  we  have  drawn  each  one  of  the  planes 
represented  in  equation  3.5  by  the  method  proposed  in 
Section  2.   The  intersection  of  the  two  planes  in  the  ct   /3 
plane  for  y   =  0  is  the  point  H  and  the  intersection  in  the 
ol   j9  plane  through  y  =   y.  ,    is  the  point  F.   Therefore,  the 


31 


^fborC-r) 


/A/TEBSECriOU  LfAJE  (PROTttXlONJ 
OKI    Qfft  PLANED 


f*or(-/3) 


Fig.  4.2  - 


Two  dimensional  representation  of  the  inter- 
section of  the  two  planes  shown  in  Fig.  4.1, 
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projection  line  is  HF.   Any  other  point  along  this  line  can 
be  determined  by:   a)  picking  up  the  corresponding  value  of 
y.     (see  y_)  and  projecting  this  value  onto  the  01   axis 
through  the  lines  on  the  Oi   y  plane  (ED  and  BA  in  this  case)  . 
These  projections  will  determine  the  points  P  and  Q  on  the 
ot   axis.   By  drawing  a  line  parallel  to  BC  through  P  and  a 
line  parallel  to  EJ  through  Q  we  determine  the  point 
y  =  y_  =  3y, ,  along  the  line  HF,  from  which  we  can  obtain 
the  corresponding  values  of  &   and  j3.   b)  Linear  interpola- 
tion can  be  used  also,  by  laying  off  the  distance  HF  along 
the  projection  of  the  line  for  equal  increments,  Ay  =  y.  . 
The  construction  method  was  used  only  as  a  graphical  help 
to  understand  the  linear  characteristics  and  the  possibili- 
ty of  linear  interpolation.   Everything  can  be  done  analy- 
tically.  The  equation  of  the  projection  line  can  be  obtain- 
ed by  using  the  formula  for  a  straight  line,  given  two 
points  H  and  F.   An  example  follows  that  will  illustrate 
the  method. 

Example  4.1:  Assume  the  system  shown  in  figure  4.3  is 
to  have  a  dominant  pair  of  poles  with  C  =  0.5  and  to  =  2.0. 
The  characteristic  equation  is 

(S  +  0.5) (S  +  1) (S  +  5) (S  +  10)+  K  K  S2  +  K  K.S  +  K  K, 

a         u       J. 

4.6 

By  rearranging  terms  we  get: 

S4  +  16. 5S3  +  (K  Ka+  73)S2  + (K  Kfc+  82.5)S  +(K  K±+    25)=  0 

4.7 
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Fig.  4.3  - 


System  with  three  variable  parameters, 
Example  4.1. 
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or 

S4  +  16.5S3+(a  +  73)S2+(£  +  82.5)S  +  (y  +  25)  =  0 

where   a  =  K  K     £  =  K  K.    y  =  K  K, 

a.  l.  1. 

by  setting   C  =  0.5   and   to  =  2.0  we  get  the  equations: 

40i  -   y  +   135  =  0 

4.8 
and  4a  -2)3  +  no  =  0 

by  making  y   =  0  we  get     cm  =   -33.75 

4.9 
and      £h  =  -12.0 

For   j3  =  0   we  get     «F  =  -27.75  4.10 

and     y   =  24.0 

We  can  now  draw  the  straight  line  on  the  &   j9  plane  as 
shown  in  figure  4.4. 

By  inspection  of  equation  4.7  we  can  see  that  one  of  the 
roots  will  go  in  the  right-half  side  of  the  s-plane  if  y 
is  negative  and  of  absolute  value  greater  than  25;  there- 
fore, for  the  system  to  be  stable  y  must  be  y  >  -2  5.  By 
setting  y  =  0  the  following  roots  are  obtained  for  equa- 
tion 4.6  : 

+ 


S 


,  2  =  -1.0  -  1.73j  (desired  value  of  roots) 


S3    =  -14.055 


S4    =  -0.446 


The  corresponding  values  of  &   and  £  are  -33.75  and  -12.0 
respectively.   Figure  4.5  shows  the  location  of  the  roots 
for  several  values  of  y.      Two  important  features  are 
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'4V_    START1M6    POINT  WITH  <K  Ab  FKEE    pAKAMeTEC. 
TOR   ROOT  lOCOS  Method  (section!    3.") 
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Fig.    4.4    -   Straight    line    representation   of    0^   line 
projection    on    the   Qlfi  plane    -   Example   4.1 


36 


noticeable  on  this  graph:   a)   Due  to  the  relationships 
established  in  equations  4„8,  the  desired  pair  of  complex- 
conjugate  roots  remains  in  the  same  place  on  the  s-plane 
for  all  values  of  y .   b)   The  remaining  roots  of  the  poly- 
nomial travel  through  the  s-plane  in  the  same  fashion  as 
on  the  root-locus.   This  is  essentially  a  method  for  deter- 
mining the  loci  of  the  roots  when  there  is  only  one  vari- 
able parameter  in  the  system.   This  fact  is  justified 
because  two  degrees  of  freedom  have  been  used  in  fixing 
two  of  the  roots  of  the  polynomial.   Therefore  it  is  pos- 
sible to  obtain  a  root-locus  for  the  remaining  roots  of 
the  polynomial.   The  method  is  outlined  in  Section  5   and 
proves  to  be  useful  in  design  for  systems  with  several 
variable  parameters.   Table  4.1  gives  the  values  of  the 
remaining  roots  and  the  corresponding  values  of  OL,    0  and  y 
for  the  characteristic  equation  of  Example  4.1. 
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Fig.  4.5  -  Locus  of  remaining  poles  in  characteristic 
equation  of  system  in  Example  4.1. 
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No. 

y 

S3 

S4 

0L 

$ 

a 

0 

-14.055 

-0,446 

-33.75 

-12.0 

b 

39 

-13.29 

-1.21 

-24 

7.0 

c 

103 

-11.77 

-2.73 

-  8 

39 

d 

231 

-7.25+3.37J 

-7.25-3.37j 

24 

103 

e 

295 

-7.25+5. 22j 

-7.25-5.22j 

40 

135 

f 

495 

-7.25+8.79J 

-7.25-8.79j 

90 

235 

g 

625 

-7. 25+10. 3j 

-7. 25-10. 3j 

122.5 

300 

h 

1007 

-7. 25+45. 3j 

-7. 25-45. 3j 

218 

491 

Table  4.1  -   Corresponding  values  of  S-  S.,    &   and  /3  for 
system  of  Example  4.1. 
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The  straight-line  property  of  the  constant  to  lines  can 

also  be  used  for  interpolation.   In  order  to  illustrate 

3 
this  let  us  assume  that  computer  results   yields  curves 

A(C  =  a,  y   ss  y   )    and  B(C  =  a,  y  =   y2)    on  figure  4.6;  but 
because  of  some  design  considerations  the  range  of  admis- 
sible values  for  Oi   and  j3  is  the  region  OPQR  on  the  QL   j3 
plane.   It  is  possible  then  to  get  a  fairly  accurate 
picture  of  the  shape  of  the  curves  in  that  region  by  draw- 
ing straight  lines  between  points  on  curves  A  and  B  having 
corresponding  values  of  to  ,  and  by  linear  interpolation  on 
these  lines  to  obtain  curves  for  intermediate  values  of  y 
between  y,  and  y_.   Curve  C  corresponds  to  y  =    (y, +  y?)/2 
and  was  obtained  by  dividing  the  straight  lines  into  two 
equal  segments.   As  it  can  be  seen  the  loop  on  curve  A 
tends  to  disappear  as  we  approach  curve  B. 
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Fig.  4.6  -  Use  of  linear  interpolation  in  constant  C 

surface  using  constant  U)  lines. 

n 
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5.   MODIFIED  ROOT  LOCUS 

There  are  many  instances  where  the  performance  speci- 
fications for  a  system  require  a  pair  of  complex  conjugate 
roots  for  the  system.   The  parameter-space  techniques  tell 
us  what  values  of  the  parameters  a,    j3  and  y   are  suitable  for 
the  system  to  obtain  this  pair  of  complex  conjugate  roots. 
Unfortunately  these  techniques  do  not  disclose  any  informa- 
tion about  the  rest  of  the  roots  of  the  characteristic 
equation  for  the  given  values  of  the  parameters  unless  we 
solve  the  polynomial,  and  if  dominance  or  stability  condi- 
tions are  not  satisfied,  we  will  have  to  go  back  and  make 
another  choice  of  values  of  Oi   fi  y.      After  several  trials  it 
is  possible  to  establish  the  trend  of  the  roots.   However, 
this  is  a  difficult  task  even  if  we  have  a  computer  program 
available   because  in  general  it  is  very  hard  to  get  a 
clear  picture,  and  in  order  to  get  detailed  results  we 
would  need  a  very  complete  three-dimensional  picture  of  the 
parameter  space.   Because  of  this  inconvenience  we  may  use 
a  method  that  gives  a  picture  of  the  system  at  one  glance 
and  tells  us  something  about  the  remaining  roots  of  the 
polynomial  which  represents  the  characteristic  equation. 
The  name  assigned  to  this  method  is  the  "modified  root  lo- 
cus", because  the  resulting  locus  follows  the  same  tech- 
niques applied  to  the  root  locus,  except  that  the  desired 
pair  of  roots  remains  in  the  same  place  on  the  s-plane 
while  the  variable  parameters  a,    /9  and  y   are  changed. 

Let  us  assume  we  have  a  system  with  characteristic 
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equation  G(S,  OL,    j3,  y)  ,    and  a  pair  of  complex  conjugate 

roots  with  given  values  of  C  and  oo  is  desired.   It  is 
^  n 

possible  to  reduce  the  order  of  the  characteristic  equa- 
tion by  making  the  division: 

g(s,  OL.    fl.  y)  K   . 

2   „-        2  * 

S  +  2£to  +  60 
n    n 

The   quotient  of    this   division  ifill  be   a  polynomial   of 

order   n-2 : 

F(S,    OL,    ft    y)  5.2 

and  the  remainder  will  be 

F1(0L,    P,    y)S    +   F2(0t,     0,  y)  5.3 

The  division  has  to  be  an  exact  one,  in  order  for 

2  2 

S   +  2Cu  +  CO   to  satisfy  the  equation  G(S,  a,    p,    y)    =0. 

Therefore,  the  remainder  (Expression  5.3)  has  to  be  equal 

to  zero,  and  for  this  condition  to  be  satisfied,  the  values 

of  F,(a,  0,  y)    and  F2(a,  j9,  y)    have  to  go  independently 

to  zero.   This  can  be  readily  seen  because  for  any  arbitrary 

value  of  S  equation  5.3  has  to  go  to  zero.   Then: 


F1(CL,     0,    y)     =    0 
and       F2^Ci>    |8'  y*  ~  ° 


5.4 


This  set  of  equations  is,  in  some  cases,  the  very  same  set 
of  equations  represented  in  equation  3. 5, or  a  linear  combi- 
nation of  them.   Therefore,   they  represent  a  line  on  the 
three-dimensional  space  OL   j9  y.   We  can  apply  the  same  tech- 
niques used  in  Section  4  to  obtain  a  graphic  or  mathematic 
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representation  and  obtain  the  corresponding  set  of  values 
of  OL,    P   and  y   that  will  yield  a  pair  of  complex  roots 
having  values  C  and  it)   .   But  what  about  the  remaining  roots 
of  the  polynomial?   Will  they  yield  the  desired  dominance 
conditions?   Will  the  system  be  stable?   What  are  the  re- 
quired values  of  ol,    j9  and  y?   Will  these  values  be  phy- 
sically attainable?   The  answer  to  these  questions  lies  in 
the  quotient  of  the  division  represented  by  Expression  5.2. 
The  fact  that  two  of  the  roots  have  been  fixed  means  that 
we  have  implicitly  used  two  degrees  of  freedom.   Therefore, 
there  is  left  only  one  degree  of  freedom  which  is  the 
essence  of  the  root-locus  method.   Its  techniques  can  be 
applied  using  the  following  procedure: 

From  equations  5.4,  we  can  get  any  two  parameters  as 
a  function  of  the  third  parameter.   Let  us  assume  we  are  to 
solve  OL   and  j3  as  functions  of  y;   then  equations  5.4  yield: 

a  =  G1(y) 

and       0  =  G2(r) , 

or,  solving  for  OL   and  y: 

OL    =    g/O) 

x  5.5 


and       y  =  G^iP)  , 


or,  solving  for  0  and  y: 

0  =  G111(a) 
and       y  =  G^1  (o>)  . 
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Now  by  substituting  the  values  of  say,  a.   and  j9  in  Expres- 
sion 5.2  (the  quotient)  we  obtain 

F(S,  y) 
or  if  we  use  &   and  y ': 

F(S,  P)  5.6 

or  if  we  use  j9  and  y: 

F(s,  a) 

Now  we  are  ready  to  apply  the  root-locus  techniques  to  the 
polynomial  F(S).  Since  this  polynomial  has  also  to  satis- 
fy the  characteristic  equation  we  can  say  that: 

F(S,  y  )    =   0  5.7 

by  rearranging  terms  of  this  polynomial  we  get 

y   N(S)  _   , 
D(S)   ~    ~    ±  5*8 

The  root  locus  for  y   will  be  the  locus  of  the  remaining 
roots  of  the  characteristic  equation.   From  this  root  locus 
it  is  possible  to  establish  stability  and  dominance  condi- 
tions.  Choose  the  more  convenient  value  of  y,      and  by  sub- 
stituting it  into  equation  5.5  get  the  corresponding  values 
of  Oi   and  j9;  or  by  a  simple  look  at  the  straight-line  method 
established  in  Section  4,  obtain  the  physically  attainable 
values  of  Oi   and  fi.      This  modified  root  locus  will  also 
allow  the  designer  to  select  the  best  choice  of  the  remain- 
ing roots  of  the  polynomial  that  will  help  to  attain  other 
performance  specifications,  such  as  bandwidth  and  gain. 
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The  method  described  above  for  three  parameters  can 
be  generalized  for  dynamic  systems  with  k  variable  para- 
meters.  In  this  case  the  designer  has  (k-1)  choices  of 
the  roots  and  there  will  be  one  "free  parameter"  that  can 
be  varied  at  will  for  application  of  the  root-locus  tech- 
niques.  The  characteristic  equation  will  be  divided: 

g(s,  a.  fi.   y k) 5>9 

Sk_1  +  b,  0Sk"2+ b 

k-2  o 

where  the  denominator  of  equation  5.9  corresponds  to  the 
product  of  the  chosen  (k-1)  roots. 

The  quotient  of  this  division  will  be  a  polynomial  of 
order  n  -  k  +  1 

F(S,  a,  j9,  y k)  5.10 

and  the  remainder  will  be 

F1(a,/3,y, k)sk"2  + Fk_2(a'^ k>s  +  Fk_1(a,^y...) 

5.11 
Again,  if  the  division  is  going  to  be  satisfied  with  the 
remainder  equal  to  zero  a  set  of  k  -  1   simultaneous  equa- 
tions are  obtained. 

Ft  (a,  j8,  y{  .  .s.  k)  =  0 


F2(a,  j3,  y   ...  k)  =  0 


Fk-2(tt'  i8'  y  ••'  k)  =  ° 


Fk_1(af  £,    y   ...  k)  =  0 


46 


5.12 


From  equations  5.12,  it  is  possible  to  solve  for  (k-1)  para- 

+--u 

meters  as  functions  of  a  k   parameter  and  by  substitution 
in  Expression  5.10  we  get: 

F(S,  free  parameter)  =0  5.13 

The  application  of  root  locus  techniques  to  equation  5.13 
will  produce  the  locus  of  the  remaining  roots  of  the  cha- 
racteristic equation.  These  techniques  are  useful  in  va- 
rious cases,  depending  on  the  applications: 

a) .   for  two  variable  parameters  it  is  possible  to  set 
one  of  the  roots  of  the  polynomial  at  a  large  negative  value 
and  obtain  a  root  locus  for  the  remaining  roots  of  the  po- 
lynomial. 

b) .   for  three  parameters  the  choice  of  the  dominant 
pair  of  complex  conjugate  roots  can  be  made,  allowing  de- 
termination of  the  remaining  roots  of  the  polynomial.   It 
is  possible  to  give  two  of  the  real  roots  large  negative 
values  and  determine  the  remaining  roots  from  the  root  locus, 
choose  the  ones  that  help  us  to  get  as  close  as  possible  to 
the  desired  performance  specifications.   Similar  considera- 
tions can  be  made  for  k  variable  parameters. 

The  calculations  involved  are  simple  and  provide  an 
overall  picture  of  the  system  and  its  possibilities.  In 
many  cases  the  quotient  is  a  function  of  only  one  of  the 
parameters;  therefore,  the  equation  for  the  root  locus  is 
available  without  the  need  of  substitutions  (see  Problem 
5.2).  With  a  computer  program  for  root  locus  the  job  is 
less  laborious. 
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Several  examples  are  available  in  order  to  clarify  the 
procedure  and  the  techniques.   Each  one  of  the  cases  out- 
lined above  is  studied. 
Example  5.1: 


OO 


KaS   +  KtS 


System  with  three  variable  parameters  -  Example  5.1. 


Fig.  5.1 
The  system  on  Figure  5.1  is  to  be  compensated  using  tacho- 
meter feedback  and  accelerometer  feedback.   We  wish  to  in- 
sure dominance  of  a  pair  of  complex  conjugate  roots  with 

C  =  0.5  and  ic     =2.0.   We  wish  also  to  know  what  are  the 
n 

permissible  values  of  the  remaining  roots  of  the  character- 
istic equations  and,  if  it  is  possible,  to  obtain  dominance 
conditions  the  required  values  of  p,  k  and  k, .  The  trans- 
fer function  for  the  system  on  Figure  5.1  is: 

100 


F(S)  = 


S5+(Of50)S4+  (SOa+SSOS^  (830CH-100i9+6250)S2+  (6250a+100rfl500)S+(L50000tf2  5 ) 

5.9 
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where   p  =  a,   K   =  0,   and  K.  =  y 

a  t 

2 

By  dividing  by   S   +  2S  +  4  we  get  the  quotient: 

S3  +  (a+48)S2  +  (48a  +  730)S  +  (730a  +  1000  +  4598) 

5.10 
and  the  remainder 

(4598a+100y  -  200^+2884)  S  +  12080a  -  400/9  -  18292 

5.11 
In  order  for  the  remainder  to  be  zero  for  arbitrary  values 
of  S,  it  is  necessary  that 

4598a  +  100y  -  2000  +  2884  =  0 

5.12 
and   1208a  -  4000  -   18292      =  0 

So  far,  two  degrees  of  freedom  have  been  used;  therefore 
we  are  left  with  only  one  of  the  parameters  variable  (free 
parameter) .   The  other  two  are  dependent  functions  of  this 
variable  parameter  (equations  5.12). 

It  is  easiest  to  have  a  as  the  variable  parameter, 
because  most  of  the  coefficients  in  the  quotient,  Expres- 
sion 5.11,  are  functions  of  a,  and  also  because  0  is 
easily  expressed  as  a  function  of  a,  in  the  second  of 
equations  5.12.   From  these  considerations: 

a        3020a  -  4573  _  , 0 

^  =  Too 5-13 

By  substituting  the  value  of  0  in  Expression  5.10  and 
making  it  equal  to  zero,  we  get: 

S3  +(a+48)S2  +  (48+730)S  +  (3750a  +25)  =0         5.14 
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Now,  by  applying  the  root  locus  technique, 

a(S2+48S+3750)       , 
S3+48S2+730S+25 

fl(S  +  24  -  56.3381)  (S  +  24  +  56.3381) 

°r  (S+23.847-J12. 377) (S+23.847+J12. 377) (S+0. 0343)  '    L 

The  root  locus  for  the  variable  parameter  &   is  shown  in 
figure  5.2.   For  any  point  on  the  given  root  locus  it  is 
possible  to  determine  the  required  value  of  01   and  the  cor- 
responding values  of  j3  and  y   are  determined  from  the  rela- 
tions expressed  in  equation  5.12.   The  graphical  method 
may  also  be  used.   The  root  locus  in  figure  5.2  shows  all 
possible  values  of  the  remaining  three  roots  of  the  poly- 
nomial.  For  low  values  of  CL,    the  real  root  will  be  very 
close  to  the  pair  of  complex  poles  that  are  desired  as  do- 
minant.  For  high  values  of  CL,    the  real  root  will  go  far 
in  the  left  of  the  s-plane  and  the  other  pair  of  complex 
roots  will  be  very  close  to  the  generating  zero  of  the 
root  locus.   As  can  be  seen,  the  system  allows  in  this  case 
an  improvement  in  the  bandwidth  by  choosing  the  value  of  CL 
that  will  make  the  real  root  lie  in  the  far  left  of  the 
s-plane  and  selection  of  the  complex  roots  to  insure  a 
dominance  condition  (real  part  of  the  complex  roots  about 
ten  times  the  real  part  of  the  dominant  roots) .   Other 
considerations  can  be  made  such  as  what  values  of  CL   will 
yield  permissible  values  of  K2  and  K, ,  without  saturating 
the  amplifiers  that  will  provide  these  gains,  etc. 
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Fig.  5.2  -  Modified  root  locus  for  remaining  roots  of 

characteristic  equation  for  system  in  Fig.  5.1 
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Example  5.2; 

A  system  has  the  following  characteristic  equation: 

G(S)  =  S5  +  16S4  +  40S3  +  (a+10)S2  +  (j3+50)S  +  (y+25) 
It  is  desired  to  have  a  pair  of  complex  conjugate  roots 
with  C  =  0.5  and  Ui     =2.0.   By  making  the  division: 

S5  +  16S4  +  40S3  +  (a+10)S2  +  (j9+50)S  +  y   +  25      5.15 

S3  +  2S  +  4 

we  obtain  the  quotient: 

S3  +  14S2  +  8S  +  (a-62)  5.16 

and  the  remainder 

(-2a+0+142)S  +  (-4a+y+273)  =0  5.17 

By  making  the  remainder  equal  to  zero  we  get: 


-2a  +    j9  +  142  =  0 
and 

-4a  +  y   +  273  =  0 


5.18 


As  is  easy  to  see  in  this  case  there  is  no  need  for  substi- 
tution, since  the  quotient  is  a  function  only  of  a.   However, 
the  values  of  ]8  and  y   will  depend  on  the  value  chosen  for  a, 
from  the  relationships  expressed  in  equations  5.18. 
Rearranging  Expression  5.16   for  root  locus  application: 

S3  +  14S2  +  8S  +  (a-62)  =  0  5.19 


or 


S3  +  14S2  +  8S  -  62 


or 

a 


5.20 


(S  -  1.74646)  (S  +  2.7266)  (S  +  13.01977) 


=  -  1       5.21 
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The  root  locus  for  equation  5.21  is  shown  in  figure  5.3. 
As  can  be  seen  it  is  not  possible  to  obtain  dominance  con- 
ditions.  The  system  will  become  unstable  for  very  large 
or  very  small  values  of  a.   In  this  case,  with  very  few 
steps  it  was  possible  to  obtain  a  clear  picture  of  the 
system.   It  is  obvious  that  the  root  locus  must  have  the 
same  shape  whatever  parameter  is  chosen.   We  can  check 
this  fact  by  choosing  say,  y  as  our  variable  parameter. 
From  equations  5.18  the  expression  for  &   is: 

«  =  y  +    273  5.22 

By  substitution  of  the  expression  for  ot   in  the  quotient 
(Equation  5.19)  we  get 

4S3  +  56S2  +  32S  +  y  +  273  -  248  =  0 
or 

1 =  -  i 

4(S3  +  14S2  +  8S  +  6.25 

If  K  =  y/4,  we  get: 

K 

(S3  +  14S2  +  8S  +  6.25) 

or 

K 


=  -  1 


(S  +  13.439) (S  +  0.28  +  jO.622) 


=  -  1  5.23 


The  root  locus  for  equation  5.23  is  shown  in  figure  5.4  and 
it  is  exactly  the  same  as  the  one  on  figure  5.3.   The  part 
which  is  not  present  is  that  corresponding  to  negative 
values  of  y.   Had  we  drawn  the  root  locus  for  negative 
values  of  y  this  part  would  appear.   The  points  where  the 
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generating  poles  for  the  root  locus  of  a  appear  are  the 
points  for  which  the  value  of  y   make  a  =   0.   Similarly 
the  points  where  the  generating  poles  for  the  root  locus 
of  y   appear  correspond  to  the  value  of  a  that  makes  y  =   0, 
We  can  make  a  quick  check.   Let  us  determine  the  value  of 
a  at  point  A  on  figure  5.3.   This  value  yields  a  =   68.25, 
and  by  substituting  in  equation  5.22  we  find  that  this  is 
exactly  the  value  of  a  that  makes  y   =  0. 
Example  5.3: 

Let  us  consider  again  the  example  4.1  in  which  the  charac- 
teristic equation  is 

S4  +  16. 5S3  +  (Q!+73)S2  +  (j3  +82.5)S  +  (y+25)  =  0 

by  using  the  same  desired  values  of  C  =  0.5  and  06     =2.0 

2 
and  dividing  by  S  +2S+4,  the  quotient  is 

S2  +  14. 5S  +  (a+40) 
and  the  remainder  is 

(-2a+j9-55.5)S  +  (-40HO/-135) 
Setting  the  remainder  equal  to  zero,  we  have 

-4a  +  y   -135  =  0 

and 

-2a  +  j8  -  55.5  =  0 


5.24 


By  comparing  equations  5.24  with  equations  4.8,  in  example 

4.1,  it  is  easy  to  see  that  these  two  sets  of  equations  are 

exactly  the  same;   therefore,  equations  5.24  represent  the 

straight  line  that  defines  a  constant  a?  on  the  constant  £ 
^  n 
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surfaces.   By  making  a  root  locus  set  up  with  a  as  a  free 
parameter  we  get: 

« =  _  ! 

(S2+14.5+40) 

5.25 

=  -  1 


(S+10.995) (S+3.705) 


The  root  locus  is  shown  as  the  slashed  lines  in  figure  5.5. 
Poles  marked  D  are  the  generating  poles  using  a  as  a  free 
parameter.   Had  y   been  chosen  as  the  free  parameter  the 
generating  poles  would  be  the  ones  marked  as  B  on  figure 
5.5,  and  the  root  locus  using  y   as  the  free  parameter  would 
be  the  part  shown  as  heavy  lines,  plus  the  slashed  lines. 
For  roots  at  points  between  B  and  D  the  value  of  &   would  be 
negative.   This  fact  can  be  checked  in  figure  4.4  where  the 
corresponding  points  B  and  D  are  marked.   The  point  of  im- 
portance in  this  case  is  that  since  root  locus  techniques 
consider  only  positive  values  of  the  variable  parameter, 
for  each  parameter  chosen  as  free  parameters,  we  will  get 
only  part  of  the  overall  root  locus.   Therefore,  if  the 
results  with  a  positive  value  of  the  free  parameter  are  not 
satisfactory,  it  is  good  practice  to  make  a  quick  check  for 
negative  values  of  the  free  parameter.   This  can  be  done 
easily  by  starting  at  the  generating  poles  and  going  in  the 
opposite  direction.   For  this  specific  example,  we  can  see 
that  for  negative  values  of  a,  the  root  locus  will  go  from 
D  towards  B^see  Figure  5.3),  and  eventually  one  of  the  roots 
will  cross  the  jcc   axis  and  the  system  will  become  unstable. 
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Fig.  5.5  -  Modified  root  locus  for  Example  5.3 
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Example  5.4-  Example  with  parameters  product. 
There  are  some  cases  where  the  characteristic  equation  con- 
tains products  of  parameters.   This  technique  is  applicable 
only  if  the  quotient  of  the  division  is  a  function  of  only 
one  of  the  parameters.   Let  us  assume  a  system  with  charac- 
teristic equation: 

S4  +  (a+15)S3  +  (20a+200)S2  +  (y+600)S  +  aj9  +  6000 

5.26 
and  the  system  is  required  to  have  a  pair  of  complex  con- 
jugate roots  with  C  =  0.5  and  U)  =10.   From  the  division 

J    3  n 

we  get  quotient 

S2  +  (a+5)S  +  (10a+50)  5.27 

and  the  remainder 

(-200a+y-400)S  +  (aj9-1000a+1000)  =  0  5.28 

The  root  locus  technique  can  be  applied   to  the  Expression 
5.27,  thus, 


or 


fl(S+10) 
(S2+5S+50) 


a(S+10) 


=  -  1 


=  -  1 


(S+2.5-6.644J) (S+2.5+6.644J) 
and  the  root  locus  can  be  seen  in  figure  5.6 
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DtSiG-eD    PAiu  OF  COMPLEX.  CON706ATE  POLES 
(BV  COINCIOEWCE  ON  THE  COOT  LOCOS  FOC  THE  PEMAi^iMe   BOOTS*) 


Fig.  5.6  -  Modified  root  locus  for  Example  5.5  -  &   as  free 
parameter. 
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The  choice  for  dominance  can  be  made  as  a  pair  of  poles  at 
S  =  -20.   From  the  root  locus  we  get 

ol  =    35 

and  from  Expression  5.28: 


a        lOOOfl  -  1000 
*  =  5 


and 


0  =  34000  = 

P  35      y/j. 


y   =  200a  +  400  =  7400 


we  could  also  obtain  the  values  by  multiplication. 

G(S)  +  (S+20)2(S2+10S+100)  =  0 

5.  29 
G(S)  =  S4+50S3  +  900S2  +  8000S+40.000  =  0 

By  comparison  of  the  two  expressions  for  G(S)  (Eqs.  5.26 
and  5.29)  we  get 

a   +  15  =  50 
or  i      ol  =   25 

2Q0L   +  200  =  900 
and        y   +   600  =  8000    :  y   =  7,400 

and       dp   +   6000  =  40.000   :  P  =    34^°  =  971 

which  checks  with  the  results  above. 

This  is  a  special  case  where  the  quotient  (Eq.  5.27) 
is  a  function  of  only  one  of  the  parameters.   In  other 

cases  the  root  locus  technique  is  not  applicable  because 

2    3 
the  free  parameter  appears  in  the  quotient  as  K,  K  ,  K  , 
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etc.,  depending  on  the  complexity  of  the  parameter  products 
in  the  characteristic  equation.   However,  there  exists  the 
possibility  that  the  designer  has  enough  freedom  to  modify 
the  system  configuration  and  make  the  selection  of  the  va- 
riable parameters  in  such  a  way  that  no  products  appear, 
or  in  case  of  products,  that  the  quotient  be  a  function  of 
only  one  parameter.   All  this  can  be  achieved  by  changing 
some  feedback  paths  or  series  compensation  schemes. 
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6 .   SINGULAR  LINES  AND  SELF-ADAPTIVE  SYSTEMS 

In  the  analysis  of  systems  with  two  variable  parameters 
by  application  of  the  transformations  proposed  by  Mitrovic 
(coefficient  plane)  and  Siljak  (parameter  plane) ,  the  pre- 
sence of  systems  with  what  is  known  as  singular  lines  was 
detected.   A  summary  of  the  procedure  that  led  to  the 
detection  of  this  type  of  system  follows.   This  leads  to  a 
better  understanding  of  what  is  meant  by  singular  lines. 
From  Section  1,  the  characteristic  equation  for  a  given 
system  can  be  written  as: 

F(S)  =  a  Sn  +  a   _ Sn_1  +  •••       +  anS  +  a   =0 
n      n-1  1     o 

6.1 
Repeating  the  manipulations  of  the  above  equation  outlined 
in  Section  1,  the  following  expressions  are  obtained 
n 

I'-^Vn*  uk-l<«>  "  ° 

k=0  6.2 

n 

!(-DkakWnk  VC)  =   0 
k=0 

Now  by  defining  the  coefficients  as 

ak  -  ^k  +  fck  +  \  6-3 

and  the  substitution  of  a,  in  Equation  6.2,  this  set  of 
equations  can  be  rewritten  as  the  system  of  linear  equa- 
tions 


ttB,  +  £c   +  D,  =  0 
0iB2    +  0C2  +  D2  =  0 


6.4 
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where  B,  ,  B2,  C,  ,  C-,  D,  and  D~  are  defined  as  in  Equa- 
tions 1.15.   When  values  of  C  and  (jo     are  assigned,  and  the 
values  of  B,,  B2,    C,,  C_,  D,  and  D2  are  calculated,  each 
of  the  equations  6.4  will  represent  a  straight  line  on  the 
Oifi   plane  and  in  a  geometric  sense,  the  intersection  of 
these  two  straight  lines  will  yield  the  values  of  Oi   and  j3 
that  satisfy  the  system  of  equations  6.4,  and  insure  that 
the  characteristic  equation  will  have  a  pair  of  complex  con- 
jugate roots  with  the  assigned  values  of  C  and  to  .   But  the 
intersection  of  two  lines  in  a  plane  is  not  always  a  point. 
Two  other  cases  are  possible:   a)   the  two  lines  are  parallel 
in  the  &P   plane  and  in  this  case  no  pair  of  a,  0  values  will 
satisfy  the  system  of  equations  6.4  (see  Fig.  6.1b),  or 
b)   the  two  lines  are  the  same  in  which  case  they  intercept 
on  an  infinite  number  of  points  and  any  point  along  this 
line  will  yield  a  pair  of  complex  conjugate  roots  with  va- 
lues C   and  0?  .   The  three  different  cases  are  shown  in  Fi- 
n 

gure  6.1.   Figure  6.1c  shows  the  case  of  coincidence  of  the 

two  lines;   this  is  known  as  a  singular  line  case  and  its 

4 
name  is  derived  from  the  matrix  algebra  manipulation   of 

the  system  of  equations  6.4.   The  systems  with  singular  line 
characteristics  find  application  in  self-adaptive  systems 
where  changes  in  the  dynamic  behavior  depends  on  the  ability 
of  the  system  to  cope  with  variations  of  the  parameters. 
The  identification  of  the  values  of  the  parameters  and  its 
changes  in  these  self-adaptive  systems  is  done  by  observa- 
tion of  the  output  of  the  system.   Identification  is  a  pro- 
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Fig.  6.1a 


£ 


Fig.  6.1b 


Fig.  6.1c 

Different  cases  of  mapping  of  a  point  on  the  s -plane 
in  the  otfi   plane. 
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blem  beyond  the  scope  of  this  work,  but  literature  on  the 
subject  is  available  (see  reference  5  and  6) .   Once  the 
change  on  the  parameters  is  identified  the  self-adaptive 
systems  use  sensitivity  models  which  are  fairly  good  for 
very  small  changes  in  the  parameters.   This  is  because  of 
the  approximations  made  in  the  selection  of  the  sensitivity 
models.   The  idea  of  using  systems  with  singular  lines  in  a 
self-sdaptive  system  came  up  because  it  yields  a  precise 
mathematical  expression  (equation  of  a  straight  line)  that 
relates  the  changing  parameters.   This  expression  can  be 
simulated  by  analog  or  digital  means,  so  the  problem  of 
steadyness  of  the  dynamic  behavior  of  the  system  is  solved 
in  part.   A  sensitivity  study  on  the  remaining  roots  of 
the  polynomial  will  tell  if  the  desired  conditions  are 
going  to  be  jeopardized  by  the  given  changes  in  the  para- 
meters and  what  tolerances  can  be  allowed  in  the  change  of 
the  parameters.   The  attractiveness  of  the  use  of  systems 
with  singular  lines  in  self-adaptive  systems  lies  in  the 
fact  that  it  is  only  necessary  that  one  of  the  parameters 
be  adjusted  for  changes  detected  in  the  other  parameter 
and  vice-versa.   As  was  said  before,  the  expression  regula- 
ting these  adjustments  is  the  equation  of  a  straight  line, 
and  if  sensitivity  studies  are  feasible,  the  tolerance  in 
the  changes  of  the  parameters  can  be  determined. 

So  far,  the  detection  of  systems  with  singular  line 
characteristics  is  a  difficult  job,  and  even  in  the  case 

where  singular  lines  are  obtained  the  values  of  C  and  0) 
^  n 
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which  generate  these  singular  lines  might  not  be  satisfac- 
tory for  the  desired  performance  specifications.   The  sub- 
ject is  relatively  new  and  research  is  going  on  on  the 
matter.   Therefore,  the  possibility  of  an  easy  method  of 
detection  and  set  up  of  singular  lines  might  appear  with 

time.   There  are  computer  programs  available  for  detection 

4 
of  singular  lines  . 

The  results  obtained  in  Section  4  and  5  of  this  work 
led  to  the  thought  that  singular  lines  on  systems  with  two 
variable  parameters  can  be  considered  as  special  cases  of 
a  system  with  three  variable  parameters  in  which  one  of  the 
parameters  is  independent  of  the  other  two  and  has  a  fixed 
value.   This  value  is  obtained  from  one  of  the  planes  re- 
presented by  equations  3.5.  It  is  clear  that  the  above  men- 
tioned plane  is  parallel  to  the  one  generated  by  the  axes 
represented  by  the  other  two  parameters. 

Figures  6.2  and  6.3  will  help  to  clarify  the  concept 
above  stated.   Figure  6.2  shows  the  general  case  where  C 
and  a:  have  been  given  and  the  point  corresponding  to  the 
given  value  of  S  maps  in  the  three-dimensional  space  as 
the  line  AB  which  is  the  locus  of  the  values  of  &,    j9  and  y 
that  will  produce  the  desired  pair  of  roots  with  parameters 
for  S  :  C  and  to  on  the  s-plane.   Line  AC  is  the  projection 
of  line  AB  onto  the  ajS  plane.   Figure  6.3  shows  the  special 
case  where  plane  A  is  parallel  to  the  aft   plane  through  y  =  y,. 
The  intercept  of  planes  A  and  B  is  the  line  AB,  that  esta- 
blishes the  relation  between  Q.   and  £  as  a  straight  line  on 
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Fig.  6.2  - 


General  case   of  mapping  of  a  point  on  the 
s -plane  into  a  straight  line  in  the  three 
dimensional  space  oifly. 
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Fig.   6.3  -  Special  case  of  mapping  of  a  point  on  the 
s-plane  with  generation  of  singular  line. 
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the  OL   &   plane  that  passes  through  y   =  y.  .   (Plane  A  in  this 
case) .   So  much  for  the  geometric  interpretation;   the  ma- 
thematical interpretation  in  terms  of  the  division  method 
outlined  in  Section  5  is  that  from  one  of  the  two  equations 
obtained  from  the  remainder,  F,  (&,  £,y)  =  0  or  F~(a,  )9,  y)  =  o 
a  specific  value  of  one  of  the  parameters  is  obtained  and 
by  substitution  of  this  value  in  the  remaining  equation,  we 
will  get  the  equation  of  the  straight  line  relating  to  the 
other  two  parameters.   An  example  follows. 
Example  6.1; 

Let  us  assume  we  have  a  system  with  characteristic  equation 
G(S)  =  S4+(lOa+lO)S3+(4Oa+50+lOy+3O)S2+(8Oa+lOj9+252) 

+  (70a+10/3+25y+340)    6.5 

and  also  assume  that  somehow  we  have  found  that  a  pair  of 
complex  conjugate  roots  with  £  =  0.5  and  to  =2.0  will 
yield  equation  F, (a, 0, y)  =  0  from  the  set  of  equations  5.4, 

defining  a  plane  parallel  to  the  OL   $   plane.   Division  of 

2 
equation  6.5  by  S   +  2S  +  4  yields  as  quotient: 

S2    +    (10«  +   8)S    +    (20G  +   5j3  +   10y   +    10)  6.6 

and   the    remainder   is 

(-20y  +   200)S   +    (-loot   -    10/8   -   I5y  +   300)  6.7 

By  letting  the  remainder  equal  to  zero  we  get 

y  -  10  =  0  6.8a 

10a  +  10)3  +  15y  -  300  =  0  6.8b 

The  graphical  representation  of  equations  6.8  is  shown  in 
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figure  6.4.   Now,  by  continuing,  from  equation  6.8a: 

y  =   10  6.9a 

and  from  equation  6.8b,  substituting  this  value  of  y: 

a  +  0  -  15  =  0  6.9b 

Now  substituting  back  in  the  characteristic  equation  the 

value  of  y   we  get 

(^(S)    =   S4+(100H-10)S3+(40a+5j9+130)  3^<8Oa+]O0+252)S  +  (7Oa  +  10&-590) 

6.10 

.  .  2 

and  dividing  by  S   +  2S  +  4  we  get  the  quotient 

S2  +  (10«+8)S  +  (20U   +  50  +110)  6.11 

Expression  6.11  is  the  same  as  6.6,  had  the  value  of  y   been 
substituted  back  in  this  expression  directly.   The  remainder 
for  the  above  division  is : 

-10a  -  10j8  +150 
and  by  making  the  remainder  equal  to  zero : 

a   +  |9  -  15  =  0  6.12 

This  is  a  singular  line  case  in  the  ot   ft   plane,  because  by 
assigning  to  y   the  value  of  10  there  are  only  two  degrees 
of  freedom  left:  a   and  $   in  equation  6.10.   According  to 
what  was  stated  in  Section  5  and  taking  G, (S)  of  equation 
6.10  as  the  characteristic  equation  of  the  system  with  two 
variable  parameters  and  using  the  two  degrees  of  freedom  in 
setting  the  value  of  a  pair  of  complex  conjugate  roots 
S  =  -1.0  -  1.732J,  the  remainder  of  the  division  should  be 
a  set  of  two  linear  simultaneous  equations  with  a  unique 
solution  for  ot   and  £. 
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Afb 


RLAHE  100(110/5 +15^-300  =0 


Fig.  6.4  -  Singular  line  case  for  system  represented  by- 
Equation  6.5  -  Example  6.1. 
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Equation  6.12  is  the  only  requirement  for  the  remainder  to 
be  zero  and  therefore  an  infinite  number  of  sets  of  values 
of  ot   and  £  can  be  assigned.   The  root  locus  for  the  remain- 
ing roots  of  the  polynomial  is  shown  in  figure  6.5,  with  ot 
as  the  free  parameter.   This  was  obtained  by  substitution  of 

/5  =  15  -  a  (Equation  6.9b) 

into  equation  6.11;  then: 

S2  +  (100!  +  8)S  +  (15a  +  185)  =  0  6.13 

or 

10tt(S  +  1.5)    i  _  . 
(S2  +  8S+  185) 

If  10a  =  K  then 


K(S+1.5) 


(S+4-J13) (S+4+J13) 


=  -  1 


By  inspection  of  the  root  locus,  and  with  the  criterion  of 
existence  of  dominance  being  that  the  real  part  of  all  of 

the  remaining  roots  of  the  polynomial  must  be  greater  or 

2 
equal  to  10  times  the  real  part  of  S   +  2S  +  4,  it  can  be 

seen  that  the  root  locus  enters  the  desired  region  at 

point  A  and  a  section  of  it  leaves  the  region  at  point  B. 

Therefore,  for  values  of  ot   between  these  two  regions  we 

have   acceptable  dominance  conditions.   For  point  A  the 

value  of  ot   is  2.41.   By  setting  ot: 

a   =  OP  +  3A  .  3^1  =  1-80 
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REGION  OF  ADM/5S/6LE  VALUES 
WQ  QEMAININ6  KOOTS 


■-+JIO 


Fig.  6.5  -  Modified  root  locus  QL   as  free  parameter 
Example  6.1  -  Singular  line  case. 
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we  can  have  tolerances  for  tt  of  +  33%  from  the  chosen 
value.   Then  for  0  we  get: 

£  =  13.19 

with  tolerance  for  /3  of  -  4.65%  from  the  assigned  value. 

It  has  been  proven  that  the  singular  line  in  the  a   fi 
plane  can  be  taken  as  a  special  case  of  the  three-parame- 
ter case.   Unfortunately,  the  detection  of  planes  parallel 
to  either  of  the  planes  containing  the  coordinate  axes  is 
a  job  as  difficult  as  the  detection  of  singular  line  in  the 
at   j9  case.   From  the  analytic  point  of  view  it  will  be  pos- 
sible to  modify  the  coefficients  of  the  different  powers  of 
S  in  such  a  way  that  a  plane  originally  at  an  angle  with 
the  planes  containing  the  coordinate  axes  becomes  parallel 
to  one  of  them.   But  from  the  practical  point  of  view  it  is 
difficult  to  get  this  modification  into  the  system,  because 
these  coefficients  come  from  manipulation  on  the  block  dia- 
gram and  usually  affect  several  terms  of  the  characteristic 
equation.   Another  limitation  may  be  that  the  desired  values 
are  not  physically  attainable,  due  to  the  intrinsic  compo- 
sition of  the  elements  of  the  system.   From  preceeding  sec- 
tions it  has  been  known  that  a  point  on  the  s-plane  maps 
in  a  straight  line  in  the  a  |3  y  space  and  this  straight  line 
has  a  projection  on  the  a   /3  plane.   Therefore,  it  might  be 
easier  from  the  designer's  point  of  view,  to "simulate"  a 
singular  line  case  by  using  three  variable  parameters,  and 
instead  of  one  set-up  for  adjustment  of  parameters,  use  two; 
that  is:  for  each  change  in  say,  y,  it  would  be  necessary 
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to  modify  accordingly  Oi   and  ft.      The  expressions  that  will 
determine  the  modification  procedure  are  precise  mathemati- 
cal expressions  represented  by  equations  5.4  and  they  can 
be  simulated  by  analog  or  digital  means  available.   How- 
ever, there  are  considerations  such  as  weight,  space  and 
economy  that  will  put  a  burden  on  a  design  of  this  type, 
but  in  the  event  that  they  can  be  overlooked,  we  will  gain 
the  ability  of  choosing  the  pair  of  dominant  poles  that 
will  help  in  the  achievement  of  the  performance  specifica- 
tions and  obtaining  solid  tolerance  limits  on  the  parame- 
ters.  Engineering  judgement  must  be  used  in  this  case  to 
weight  all  these  considerations. 

When  the  design  considerations  stated  above  become  a 
real  problem  and  a  singular  line  case  is  desirable,  it  is 
worth  while  to  study  the  coefficients  of  the  polynomial 
which  represents  the  characteristic  equation,  and  by  alge- 
braic manipulation,  obtain  the  singular  line  case.   The 
method  outlined  in  the  following  section  for  a  fourth  order 
polynomial  can  be  used  in  higher-order  polynomials.   The 
results  obtained  are  quite  satisfactory  and  a  clear  picture 
of  what  has  to  be  done  to  obtain  singular  lines  is  presented, 
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6.1   GENERATION  OF  SINGULAR  LINES  IN  A  FOURTH-ORDER 
POLYNOMIAL  WITH  THREE  VARIABLE  PARAMETERS. 
A  Fourth-Order  characteristic  equation  of  a  system 

with  three  variable  parameters  can  be  written  in  general 

form  as : 


S4   +    (b3a+c3j3+e3y+f3)S3   +    (b20+c2j9+e2y+f2)S2   + 


6.14 


(bna+cn  j9+e.,y+f)S   +    (b   QH-c    0+e   y+f   )    =  0 
1        1        1  v    o        o        o'      o' 

By  making    the   substitution   of 

\   =   V  +   fk  6'15 

equation   6.14   can  be  written   as: 

S4    +    (b3a+c3/3+d3)S3    +    (b2a+c23+d2)S2    +    (h^+c^+d^  S    + 

(b   a+c    0+d    )    =   0  6*16 

v    o        o        o 

It  is  clear  that  the  case  of  singular  lines  will  depend  on 
the  configuration  of  the  polynomial,  that  is,  on  the  coeffi- 
cients of  the  different  powers  of  S.   Therefore,  it  is  pos- 
sible to  determine  what  coefficients  must  be  affected  in 
order  to  obtain  singular  lines.   If  it  is  assumed  that  the 

polynomial  has  singular  lines  for  values  of  x  =  2Cco  and 

2  2 

y  =  w     ,    division  of  the  polynomial  by:  S   +  x  S  +  y  yields 

the  follwoing  results: 

Quotient  s 

2  2 

S   +  (b30H-c3j3+d3-x)S  +  (b2-xb3)a  +(c2-xc3)£  +(d2-y-xd3+x  ) 

6.17 
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The  remainder  is  j 


F1(Oi,  j3)S    +   F2{OL,$)  6.18 


where : 


2  2 

F,  (a,  0)    =    (b,-yb    -xb2+x  t>3)tt  +    (c,-yc3-xc2+x   c_)  0  + 

(d1-yd3+2xy-xd2+x2d3-x3)  6-19 


and 

F2(a, 0)    =    (bQ-yb2+xyb3)a  +    (co-yc2+xyc3) 0  + 

2  2  6,2° 

(dQ-yd2+y   +xyd3~yx    ) 

Singular  lines  can  be  generated  for  each  of  the  following 

cases  : 

Case  I :   F,  (a,  0)  =  0  for  arbitrary  values  of  Oi   and  0.  Then: 

bl  "  yb3  "  Xb2  +  x2b3  =  °  \     6.21 

2 

c,  -  yc.  -  xc?  +  x  c   =  0  6.22 

2       3 

d,  -  yd_  +  2xy  -  xd?  +  x  d_  -  x   =  0  6.23 

which  is  a  system  with  three  equations  and  two  unknowns  (x 
and  y) .   From  inspection  of  the  system  the  following  conclu- 
sions are  drawn : 

a) .   It  is  possible  to  solve  for  x  and  y  by  utilizing 
only  equations  6.21  and  6.22. 

b) .   Equation   6.23  is  a  function  only  of  d, ,    d2  and 
d..,  which  in  turn  are  functions  of  y    (equation  6.15). 
Therefore,  once  the  values  of  x  and  y  are  found,  y   can  be 
set  at  the  value  that  will  satisfy  equation  6.23  for  the 
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values  obtained  for  x  and  y.   It  is  worth  noting  that  the 
modifications  in  equation  6.23  do  not  affect  the  values  of 
x  and  y  obtained  before.   Then,  from  equations  6.21  and 
6.22: 

9-  Clb3  -  blC3 

2l,U)     =  x  =  — r r 

c2b3  -  b2c3 

6.24 

2 
9       b,  -  xb„  +  x  b_ 

to   =  y  =  r 

n  b3 

and  the  third  equation  to  be  satisfied  once  x  and  y  have 

been  found  is : 

2       3 
d.  ~  Y^-5  +  2xy  -  xd9  +xd_-x   =0        6.25 

By  substituting  d,  =  e,  y  +   f,  and  using  the  values  obtained 
for  x  and  y  we  can  solve  for  y. 

If  the  values  of  C  and  U)     are  not  satisfactory  to  the 
designer,  it  is  possible  to  know  exactly  what  coefficients 

can  be  modified  in  order  to  get  the  desired  values.   In  this 

2 
case  if  for  example  y  =  o>   is  satisfactory  but  we  wish  to 

change  x  =  2£to  ,  then  C,  ,  C2,  b,  and  C-.  can  be  modified  at 

will  without  disturbing  the  value  of  y.   From  the  above 

coefficients  the  designer  might  be  able  to  choose  the  most 

convenient  methods  for  the  acutal  modification,  such  as 

feedback  paths,  series  compensation  devices,  etc.   The 

above  considerations  are  applicable  to  cases  II  and  III 

which  follow. 
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Case  II.  F0(Oi,  P)    =  0   for  arbitrary  values  of  a.   and  j3.  Then 


b   -  yb_  +  xyb_  =  0 
o      2  3 


6.26 


Co  "  yC2  +  Xyc3  =  ° 

2 

lo  "  yd2  +  Y   +  Xyd3 


2  2 

d   -  yd0  +  y   +  xyd_  -  yx   =  0 


From  the  first  two  equations  we  get: 

b  c~  -  b«c 

~v».  ,  o  2 2  o 

2C^  =  x  =  r 7- 

n       b  c_  -  b-c 
o  3     3  o 


2  b 


6.27 


and 


n    J    b?  -  xb_ 


2  2 

d   -  yd2  +  y   +  xyd   -  yx   =  0  6.28 


Again,  by  substitution  of  d,  =  e,  y  +  f  ,  and  using  the 
values  obtained  for  x  and  y,  it  is  possible  to  solve  for  y, 
Case  III.   F,  (a,  j3)   =   F~(a,  j3)  for  arbitrary  values  of  a 

and  0.   Then: 

2 

b,  -  yb_  -  xb~  +  x  b_  =  b   -  yb0  +  xyb_ 
1-^3      2       3     o    -*  2     -*  3 

2 

c,  -  yc  _  -xc0  +xc_  =c   -  yc  ~  +  xyc _  6.29 
13      2       3     o      2     ^3 

2       3  2  2 

d,  -  yd^  +  2xy  -  xd2  +  x  d_  -  x   =  d  -yd2+  y  +  xyd_-yx 

From  the  first  two  equations  we  get 

(b2  -  b3)  (Cl  -  cQ)  +  (bl  -  bQ)  (c3  -  c2) 
n       b3(Cl  -  c2  -  eo)  +  c3(b2  -  h1   +  bQ) 

x.     -u   ■   2,     .  6.30 

~       b,  -  xb_  +  x  b_  -  b 
, ,  2        1 2 3 g 

%  =y =  b.  +  xb.  -  b. 
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and 

[1 


2       3  2  2 

d,  -  yd   +  2xy  -  xd0  +  x  d   -  x   =d-  yd0  +  y  +  xyd~-  yx 


6.31 
That  with  the  given  values  of  x  and  y  will  yield  the  value 
of  y. 

The  important  feature  of  the  results  stated  above  is 
that  the  designer  knows  exactly  what  coefficients  of  the 
characteristic  equation  affect  directly  the  values  of  C 
and  cc     and  hence,  what  parts  of  the  actual  system  must  be 
modified  in  order  to  obtain  the  desired  values. 

The  solutions  proposed  above  are  for  the  case  where 
the  "free  parameter"  is  y;   that  is:   one  of  the  planes 
considered  in  Section  6  is  parallel  to  the  &P   plane  through 
the  calculated  value  of  y.   The  quotient  (equation  6.17) 
will  be  then  a  function  of  Oi   and  j9.   From  the  remainder,  a 
can  be  expressed  as  a  function  of  j9,  or  0  as  function  of  a, 
and  root  locus  techniques  applied  in  order  to  determine  the 
location  of  the  remaining  roots  of  the  polynomial.   A  sensi- 
tivity study  can  be  made  on  the  root  locus.   For  applica- 
tions in  self-adaptive  systems  there  will  be  need  of  only 
one  set-up  for  modifications  when  changes  in  either  of  the 
parameters  occur.   The  following  example  will  illustrate 
the  technique. 
Example  6.2; 

A  system  has  the  following  characteristic  equation: 
S4  +  (10a+10)S3+  (40a+5j3+130)S2+  (80a+10j9+10y+5)  S  + 

(70a+10j9+y+200)  6.32 
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Question:   What  values  of  C  and  to  are  attainable  in  sin- 
gular line  cases  and  what  are  the  required  values  of  y? 

For  Case  I,  substituting  the  values  of  the  coefficients 
in  equations  6.24: 

x  =  2Cw  =2    or    C  =  0.5 
n 

2 
y  =   0)   =4    or   0)  =  2  6.33 

J  n  n 

y   =  24.7 

Substituting  the  value  of  y   in  equation  6.32  and  dividing 

2 

by  S   +  2S  +  4   we  get: 

Quotient: 

S2  +  (100H-8)S  +  (208+50+110)  =0  6.34 

Remainder: 

-100!  -100  -  215.3  =  0  6.35 

Now  root  locus  techniques  can  be  applied  to  equation  6.34 
after  substituting  0  as  a  function  of  0L   from  equation  6.35. 
The  expressions  for  the  quotient  and  the  remainder  can  also 
be  obtained  from  equations  6.17  and  6.20  simply  by  substi- 
tuting the  values  of  the  coefficients. 

For  Case  II,  substituting  the  values  of  the  coefficients 
in  equations  6.27: 

x  =  2C^  =  0.5    or    C  =  0.1765 

y  =  oo  2    =   0.2    or   0)  =  1.4142  6.36 

2  n  n 

y  =   46.5 
Substituting  the  value  of  the  y   in  equation  6.32  and  divid- 
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2 

mg  by  S   +  0„5S  +  2   we  get: 

Quotient: 

S2  +  (100H-9.5)S  +  (353+50+123.25)  =  0  6.37 

Remainder: 

(42.5tt  +  7.5j9  +  389.375)S  =  0 

For  Case  III,  substituting  the  values  of  the  coefficients 
in  equations  6.30: 

2£<o  =  x  =  1   or    C  =  0.5 
n 

2 
a>   =y  =  l   or   0)=1  6.38 

n    J  n 

y  =   204/9  =  22.667 

Division  yields 
Quotient: 

S2  +  (10tt+9)S  +  (30QH-5J9+120)  =  0  6.39 

Remainder: 

(40OH-5j3+102.667)S    +    (40QH-5i9+102 .  667)    =   0  6.40 
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7.   CONCLUSIONS  AND  RECOMMENDATIONS  FOR  FURTHER  WORK 

In  this  work  it  is  shown  that  the  increase  in  degrees 
of  freedom  in  a  system  really  helps  the  designer  in  the 
task  of  obtaining  from  it  the  desired  characteristics. 
Mitrovic's  approach  has  proved  to  be  useful  in  the  parame- 
ter plane  or  when  products  of  parameters  are  present  in 
the  characteristic  equation  of  a  system.   However,  when  the 
number  of  variable  parameters  is  three  or  more  it  is  diffi- 
cult to  obtain  graphical  information.   The  tabulated  data 
obtainable  from  a  computer  comes  only  in  quantized  form  and 
it  is  difficult  for  the  designer  to  use  it  in  seeking  his 
objective,  which  is  to  make  the  best  choice  of  the  variable 
parameters.   In  the  field  of  singular  lines  a  method  of  ge- 
neration of  singular  lines  for  a  fourth-order  polynomial 
was  established  and  it  is  expected  that  formulae  can  be 
obtained  for  systems  with  characteristic  equation  of  higher 

order.   With  this  method  the  values  of  C  and  0)  at  which 

n 

singular  lines  are  expected  and  the  value  at  which  a  third 
parameter  must  be  set  are  obtained.   In  addition,  the  mo- 
difications in  the  system  for  which  the  values  of  C  and  Ui 

•*  n 

are  changed  are  indicated. 
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